Feeding upon negative entropy in a thermal-equilibrium environment 



Bruno Crosignani 

Department of Applied Physics, California Institute of Technology, Pasadena, California 91125, USA and 
Dipartimento di Fisica, Universita' dellAquila, 67010 LAquila, Italy and 
Istituto Nazionale per la Fisica della Materia, Universita' di Roma "La Sapienza" , 00185 Roma, Italy 

Paolo Di Porto 

Dipartimento di Fisica, Universita' dell'Aquila, 67010 L'Aquila, Italy and 
Istituto Nazionale per la Fisica della Materia, Universita' di Roma "La Sapienza" , 00185 Roma, Italy 

Claudio Conti 

Istituto Nazionale per la Fisica della Materia, Universita' Roma Tre, 00146 Roma, Italy 

(Dated: February 9, 2008) 

The validity of the Second Law of thermodynamics, indisputable in the macroscopic world, is 
challenged at the mesoscopic level: a mesoscopic isolated system, possessing spatial dimensions of 
the order of a few microns, is capable, as shown by a straightforward kinetic analysis, to exhibit a 
perpetuum mobile behavior associated with large negative variations of the Clausius entropy of the 
system. This violation of the Second Law is expedient for devising a cyclic process through which 
an isolated system can extract energy from a surrounding thermal bath. 



PACS numbers: 

As first clearly stated by Schroedinger "a living organ- 
ism tends to approach the dangerous state of maximum 
entropy, which is death. It can only keep aloof from 
it, i.e., alive, by continually drawing from its environ- 
ment negative entropy. What an organism fed upon is 
negative entropy". [l| Is it possible for an organism im- 
mersed in a thermal bath to be insulated, and thus avoid 
thermalization, and still be able to reduce its entropy? 
Unfortunately, the above scenario contradicts one of the 
most cherished laws of physics, that is the Second Law 
of Thermodynamics, which has victoriously resisted all 
the attempts aimed at finding a particular case where it 
could be violated. However, one has to consider that the 
Second Law has been formulated in the context of macro- 
scopic physics and it is in this context that it has been 
successfully applied and verified. While the Second Law, 
being inherently statistical in its nature, cannot carry 
over to microscopic cases where the number of involved 
particles is too small, its limits of validity are not well 
understood in the region bridging the macroscopic to the 
microscopic world. In particular, the intermediate meso- 
scopic regime, which includes part of the biological realm, 
is still terra incognita and open to possible surprises. Re- 
cent years have actually witnessed a growing interest in 
the thermodynamics of small-scale non-equilibrium de- 
vices, especially in connection with the operation and ef- 
ficiency of Brownian motors. Q We believe that we have 
been able to devise a particular case in which an isolated 
object, placed in a thermal-equilibrium environment, is 
nevertheless capable of reducing its entropy or, equiva- 
lently, to extract work from the thermal bath surround- 
ing it. We have identified this object with a mesoscopic 
system, typically of the size of a few microns, consist- 
ing of two isolated cavities, filled with a common number 



of molecules, separated by a movable adiabatic imper- 
meable partition (see Fig^). Its macroscopic equivalent 
is the so-called adiabatic piston, sometimes referred to 
as "enigmatic" because the description of its dynamical 
evolution toward equilibrium is far from being trivial, the 
problem being actually undetermined from the point of 
view of elementary thermodynamics. 3] A careful analy- 
sis of the problem, which has been tackled a few years 
ago by an approach based on gas kinetic theory, shows 
that thermodynamic fluctuations cannot be neglected at 
the mesoscopic level and that they actually play a major 
role in determining the dynamic evolution of this isolated 
system. As shown in this paper, its behavior implies that 
the system fails to settle down to a well-defined equilib- 
rium state and undergoes negative entropy fluctuations 
whose ensemble average over many replicas of the sys- 
tem can be quite large. Thus, on the average, the system 
is able to feed on negative entropy. It is also possible 
to show how the system can function as an engine and 
turn the entropy extracted from the surrounding thermal 
bath into work. These remarkable results are derived in 
the frame of a model which does not contain any ad hoc 
parameters and appears to possess some simple universal 
feature which makes it a potential candidate for support- 
ing some still unknown ubiquitous biological process. 

Our system consists of two isolated cavities A and B, 
filled with a common number of moles of the same gas, 
separated by a movable impermeable partition which is 
assumed to be initially at rest in a position corresponding 
to equal volumes and temperatures (and thus pressures) 
in the two cavities. Under these conditions, one would 
expect the wall to remain at rest, apart from microscopic 
random oscillations around the equilibrium position. 

This is obviously the case in most ordinary situations, 
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FIG. 1: (Colors online) The adiabatic piston: an insulating 
cylinder divided into two regions by a movable, frictionless 
and insulating piston. Left (thin line): evolution as a function 
of t — t/to of the rms deviation of the piston position from 
its initial value £(0) = 1/2, for fi = 0.5 and N = 3 • 10 4 (over 
1000 realizations). Right (thick line): time evolution of the 
entropy. 



and pressures in the two sections. Once the piston is re- 
leased, the system evolves in a labile quasi-equilibrium 
state, driven by the random motion of the wall induced 
by the elastic molecular collisions, while the gas pressure 
on both sides remains equal to the initial one. 0, Q We 
describe the dynamics of the system in terms of the ran- 
dom variable X(t), that is the instantaneous position of 
the piston, whose ensemble-average value (X(t)} remains 
equal to L/2 for obvious symmetry reasons. The normal- 
ized mean-square value ([X(t)~L/2] 2 )/(L/2) 2 represents 
the ensemble- aver age of the quantity (AV) 2 /(Vb/2) 2 ap- 
pearing in Eq.(l) and determines the corresponding en- 
tropy decrease undergone by the system. 

The system evolution from an initial configuration of 
macroscopic equilibrium, that is equal pressures and tem- 
peratures on both sides and piston at rest, is described 
by the stochastic equation[5j 



where the wall thermal conductivity is large enough to 
maintain a common temperature in the two cavities, so 
that any small random displacement of the wall is imme- 
diately reversed by the induced pressure change. Con- 
versely, if thermal conductivity is negligible, it will be 
shown below that the wall can actually undergo size- 
able random displacements, during which the pressures 
on the two sides remain the same. In each realiza- 
tion, the system evolves and the total entropy variation, 
from the initial state in which both cavities have the 
same volume Vq/2 and temperature To to a generic one 
(Va,T a ;V b ,T b ), is 
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where AV — Va — Vq/2, c p being the molar heat at con- 
stant pressure and n the common number of gas moles 
on each side. The ensemble average AS of this intrinsi- 
cally negative entropy change can assume large values in 
specific mesoscopic situations, as it will be proved in this 
paper. Note that in deriving Eq.QJ, we have exploited 
the relation Ta + Tb = 2T , which corresponds to ne- 
glecting the kinetic energy of the wall. In fact, under 
the adiabatic hypothesis, it can be regarded as a body 
possessing a single (translational) degree of freedom. As 
a consequence, its energy and entropy are of the order 
kTo/2 and k, respectively, and can be neglected. 

In order to corroborate the above considerations, we in- 
vestigate the so-called "adiabatic-piston problem" , 0, Q 
dealing with a system which, in the language of classical 
thermodynamics, consists of an isolated cylinder divided 
into two parts by means of a frictionless adiabatic piston, 
each section containing the same number n of moles of 
the same perfect gas (see Fig.JTJ). Initially, the piston is 
held in the fixed position corresponding to equal volumes 
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where k is Boltzmann's constant, N the common num- 
ber of molecules on each side, and /x = M/M g , M and 
M g being respectively the mass of the piston and the 
common value of the gas mass in each side. Deter- 
mining the correct expression of a(t) is a delicate task, 
since the standard Langevin approach does not in gen- 
eral carry over to nonlinear dynamical systems 0, as 
the one described by Eq.JSJ). In order to take advan- 
tage of the Langevin method, we linearize the above 
equation: a) by considering small displacements around 
the starting position, that is \(X - L/2)|/(L/2) < 1, 
and b) by approximating the square of the piston ve- 
locity (dX/dt) 2 with its mean-square velocity kTo/M 
(both hypotheses will be proved consistent a posteriori) . 
Proceeding in this way, a straightforward application of 
the dissipation-fluctuation theorem yields (a(i)a(t')) = 
[8(2m/cT/7r) 1 / 2 PS , /M 2 ]<5(t - t'), where m is the mass of 
the individual molecule, P the common pressure on the 
two sides of the piston and S its area. After introducing 
the variable x = X — L/2, the linearized form of Eq.J2J) 
reads 



a(t). (3) 



We now observe that the above equation is formally iden- 
tical to the one describing the Brownian motion of a 
harmonically-bound particle of mass M, that is 




x + fix + uj 2 x = A(t), 



(4) 



where A(t) is the Langevin acceleration, a problem 
which has been thoroughly described in the literature 
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7J. By comparing Eqs.JSJ) and (@J), we can obviously ap- 
ply the results of to our case by identifying j3 with 
8(2NkT Q /iriiML 2 ) 1 / 2 and uj 2 with 8kT Q /fiML 2 . When- 
ever (3 3> u> ( "overdamped" case) , which in our situation 
is equivalent to the obviously satisfied relation y/~N ^> 1, 
the analysis carried out in naturally highlights the 
existence of two significant time scales tth — 1/2/3 and 
tas = /3/2u 2 (t as ^> t t h)- They represent the thermaliza- 
tion time tth, i-e., the time over which the mean-square 
velocity ((dx/dt) 2 ) attains its equipartition value kTo/M, 
and the much longer time t as over which the mean- 
square displacement reaches its asymptotic value (a; 2 ) = 
kT /Mui 2 . For our system t as = (NL/w)(M/M g )/TT 1 ^ 2 
and t as /t th = IQN/n, where w — (2kT Q /m) 1 ^ 2 is the 
most probable velocity of the gas Maxwellian distribution 
function. Therefore, since N ^> 1, tth, is much smaller 
than t as , a circumstance which justifies a posteriori the 
replacement of (dX/dt) 2 in Eq.© by its average value 
kTo/M. The asymptotic value of the mean-square dis- 
placement of the piston from its central position reads 



< x > = 



kT = (i L 2 
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(5) 



so that, in the limit < 1 (that is, small piston mass 
with respect to gas mass), the above assumptions a) and 
b) allowing us to linearize Eq.@ are both satisfied. 

We can now evaluate the entropy change (AS) aver- 
aged over many realizations of our system, starting from 
the initial state corresponding to X(t — 0) — L/2. Since 
in our case ((AV) 2 )/(V /2) 2 = (x 2 )/(L/2) 2 = fi/2 < 1, 
Eq.Q) approximately yields 
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k 
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which corresponds to a large entropy decrease whenever 
[iN 3> 1. The question naturally arises: what are the 
spatial and temporal scales over which this violation of 
the second law of thermodynamics can actually occur? 
Before answering this question, the conceptual meaning 
of the result implied by Eq.© has to be clarified. To this 
aim, we note that, according to standard thermodynam- 
ics, any closed system in an equilibrium state, once an 
internal constraint is removed, eventually reaches a new 
equilibrium state characterized by a larger value of the 
entropy. The case considered in this paper represents a 
remarkable exception to this statement. More precisely, 
our closed system is made up of the two gases and the 
initial internal constraint is provided by the piston held 
in the central position. After the piston is released, our 
analysis shows that no final equilibrium state is even- 
tually reached, since the piston keeps wandering, per- 
forming random oscillations around X=L/2, with mean- 
square amplitude (x 2 ) — (L/2) 2 fi/2. 

The system behaves as a perpetuum mobile of the sec- 
ond kind. If embedded in a thermal bath at temperature 
Tq and pressure Pq, one can devise a process through 
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FIG. 2: Cyclic process through which work can be extracted 
from a thermal bath at temperature To. 



which work can be extracted from this environment. To 
this end, let us assume our system to be capable of split- 
ting into the two separate cavities A and B if a large 
fluctuation AV^ occurs (see Fig. EJ. The two cavities 
can then undergo a reversible adiabatic process which 
drives them back to the initial temperature Tq. If, at 
this point, the wall of the cavities becomes thermally 
conductive for a finite amount of time, the two cavities 
can be brought back to the initial volume Vo/2 through 
a reversible isothermal process at the bath temperature 
To. Let us now assume that the wall recovers its insulat- 
ing nature and that the two cavities are reunited again: 
the total work extracted in this cyclic process is precisely 
W = T \AS\ = k Tq Nfi, where A S is the negative en- 
tropy variation associated with A V. 

We now return to the question of the spatial and tem- 
poral scales over which our results apply. To this end, 
let us consider the specific case of a gas under standard 
conditions of temperature and pressure, for which, ex- 
pressing hereafter L in microns, N = 3 • 10 7 L 3 . Accord- 
ing to Eq.@, if we refer as an example to a biatomic 
gas, we have (AS)/k = —5 • 10 7 L 3 /i and, recalling the 
expression of i as , that is t as = NL/i/wt/tt, we obtain 
t as = 5 ■ 10 _2 /xL 4 sec (having assumed w = 4 • 10 s mi- 
crons/sec, molecular oxygen). This extremely sensitive 
dependence of t as on the linear dimension of the cylinder 
(t as (x L 4 ) appears to limit the applicability of our model 
to values of L up to a few microns. In fact, beyond this 
mesoscopic scale, t as becomes so large as to render un- 
realistic the adiabatic piston assumption over this time 
interval. As an example, for fi = 10~ 2 and L = 1 cm, 
we obtain t as = 5 ■ 10 12 sec, that is about 1000 centuries! 
Conversely, by taking L = 1 micron, we get the reason- 
able value t as = 5 • 10" 4 sec and (AS) /k = -5- 10 5 . This 
corresponds a violation of the second law in the meso- 
scopic realm. 

The above approach has allowed us to deal with the 
situation -C 1. In order to have an insight into the 
behavior of our process in the more general case fj, < 1, we 
assume Langevin's approach to be approximately valid 
also in this moderately nonlinear regime, and use the 
nonlinear Eq.Q with the same stochastic acceleration 
a(t) worked out in the linear case. After introducing 
the dimensionless units £ = X/L and r = t/t a , where 
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FIG. 3: Entropy change for r = 10 6 as a function of /i and 
N = 3 • 10 4 (over 1000 realizations). The continuous line 
represents the predicted behavior in the linear regime. 

t Q = A^/2t th , Eq.© reads 

(7) 

where the dot stands for derivative with respect to r, 
a(r) is a unitary-power white noise process and a 2 — 
(7r/2\/2)0/iV). This last equation can be numerically in- 
tegrated by adopting a second-order leap-frog algorithm 
as the one developed in Q. In particular, we are able 
to evaluate the ensemble average (over 1000 realizations 
of the piston) of the time evolution of the mean square- 
root deviation (x 2 ) of the piston position. This analysis 
is reported in Fig.lQJ for /i = 0.5 and N = 3 • 10 4 . Its 
inspection clearly shows that the piston undergoes ran- 
dom fluctuations around the central position X = L/2, 
which increase with time up to an asymptotic value of 
the order of L/10. The ensemble average of the total 
entropy change (AS)/k = N(cp/R)(ln(l - 4x 2 /L 2 )) of 
our isolated system can also be numerically evaluated. 
In Fig.JQ), (AS}/k is reported as a function of the nor- 
malized time r. 

Further numerical analysis shows (see Fig.©) that, 
keeping N fixed and varying fi, the asymptotic value of 
(AS) jk undergoes an approximately linear decrease with 
fx, up to /i = 0.6, as predicted by Eq.JJjJl. Beyond this 
value the relative entropy decrease slows down, but we 
are now in a range where Langevin's model may not be 
valid. Conversely, in the limit /i — > our approach fur- 
nishes (X- < X >) 2 /(L/2) 2 = n/2 -> and t as -> oo, 
results which are consistent with those of 

We wish to note that the celebrated Boltzmann's H- 
theorem ^(Ji which can be considered as a "proof" of 
the second law, does not apply in our case since one of 
its main hypotheses, that is the complete molecular-chaos 
assumption (Stosszahlansatz) , is not valid in our dynam- 
ics: in fact, the correlation induced by the random mo- 
tion of the piston favors a common sign of v x and v' x in 



the two-particle correlation function of the gas near the 
piston, while the single-particle correlation functions are 
independent from the sign of v x and v' x . This prevents the 
two-particle correlation distribution to factorize into the 
product of the one-particle distribution functions and, 
thus, the molecular-chaos assumption does not apply. 

Finally, the validity of our conclusions appears to be 
corroborated by suitable molecular dynamic simulations 
of the evolution of our system. These typically involve a 
considerable number of point particles, which model the 
gas inside the cylinder, separated by a frictionless piston 
against which they undergo perfect elastic collisions. In 
particular, a microscopic model consisting of N = 500 
hard disks furnishes a relaxation time r (corresponding to 
our t as ) in good qualitative and quantitative agreement 
with out results. [llj 

Besides, numerical investigations of a system consist- 
ing of a number of particles of the order of 10 3 , indi- 
cates that the difference between the temperatures on 
the two sides undergoes relevant oscillations, so that the 
system does not reach equilibrium, as predicted by our 
model. [3 

Further recent molecular dynamic simulations describ- 
ing the time evolution of the piston position around 
X = L/2 are in fairly good agreement with our results, 
both qualitatively and quantitatively. [13| 

We wish to thank Noel Corngold for many fruitful dis- 
cussions. 
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